Quantum discord and non-Markovianity of quantum dynamics 
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The problem of recognizing (non-)Markovianity of a quantum dynamics is revisited tfirough analyzing quan- 
tum correlations. We argue that instantaneously-vanishing quantum discord provides a necessary and sufficient 
condition for Markovianity of a quantum map. This is used to introduce a measure of non-Markovianity. This 
measure, however, requires demanding knowledge about the system and the environment. By using a quantum 
correlation monogamy property and an ancillary system, we propose a simplified measure with less require- 
ments. Non-Markovianity is thereby decided by quantum state tomography of the system and the ancilla. 
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I. INTRODUCTION 

Real quantum systems are inevitably open because of in- 
teraction with their ambient environment. As a result, in- 
formation, in principle, can leak into the environment, and 
may be fed back later in a different form into the system 
Jill. This scenario makes accurate description of the related 
dynamics difficult ||2J], because one often has no much in- 
formation/control about/over the environment — unless, e.g., 
the environment can be controllably engineered fl. Never- 
theless, under some specific (ideal) conditions, such as weak 
coupling with a memoryless environment and Born-Markov 
approximation, the underlying system evolution can be well 
described by a quantum Markovian dynamical equation [ill 
(see also Ref. |01 for an alternative derivation). This approx- 
imation has been proved useful in various situations, such as 
devising quantum error-correction schemes [Q . 

In this case, the dynamics of the system for any interval 
T — regardless of the start time to — is given by a completely- 
positive (CP) map £{t), having the "dynamical semigroup" 
roperty £{t)£{t') = £{t + t'), for any t,t' > 
Sy; or equivalently, the dynamical equation dg{T)/dT = 
£[p(r)] governs the system density matrix, where the time- 
independent C has the following Lindblad form: 

£[•] = -i[H, ■] + ^r™(2F„ • F;„ - {F;„F„, •}+),(!) 

m 

leading to £{t) = e"^^. Here, r^ > Vtti, H is a Hermitian 
operator, and Fm& are some operators acting on the system's 
Hilbert space j^^. Note that the case of time-dependent r^s 
is also referred to as a time-dependent Markovian evolution 
whenr,„(f) > OVtH. 

This simplified picture (even its beyond-Markovian 
version — Redfield theory), although fairly well applicable to 
various situations, often fails to capture dynamics of complex 
quantum manybody systems in condensed matter or in some 
biological complexes (which may feature quantum effects). 
In such cases, at least one of the assumptions of the Bom- 
Markov (or Redfield) theory breaks down llTol - [T2l . 

These difficulties have spurred introducing various for- 
malisms to somehow incorporate non-Markovian effects into 
open system's dynamics, e.g., through adding memory ker- 
nels 1I13I1 or even devising specific dynamical equations valid 
for non-Markovian regimes (see, e.g., Refs. llll[l4ll). 



Although deciding whether a given quantum chan- 
nel/dynamic^Js Markovian (or not) is in complex-theoretic 
sense hard flSll . there exist a number of measures of Marko- 
vianity based on the Lindblad form JaQ]. In addition, some 
experimentally measurable criteria have been introduced for 
detecting non-Markovianity of an evolution by using, e.g., 
increasing (decreasing) the distance (fidelity) of two density 
matrices in time llmilTll . or increasing the entanglement be- 
tween the system and an isolated ancilla in some instances of 
time II18I1 . These measures provide sufficient conditions for a 
channel to be non-Markovian, which leaves open search for 
more delicate measures — for a critical study and comparison 
of some of these measures, see Ref. |19']. 

On a related note, since correlations and how they are dis- 
tributed between a system and its embedding environment 
play a crucial role in the Markovianity property of the re- 
sultant system dynamics, it seems natural that tracking cor- 
relations may provide a hint for (at least partially) identi- 
fying Markovianity property of the dynamics. Along this 
line, here we employ a recently introduced quantum correla- 
tion measure — "discord" Ii20il — and its associated properties 
to propose a (non-)Markovianity measure. We illustrate the 
main idea through two examples. 



II. QUANTUM DISCORD AND MARKOVIANITY 

Quantum discord (QD) is a measure for quantumness of 
correlations in a bipartite state qse ||20M22I1 . and is defined as 

VsIqse] = S[gs] + min S [S|{nf }] - S[gsE], (2) 

where S is the von Neumann entropy, and iS[ii^|{n:?}] = 
^ PjS[gEiYi^] is the average entropy of the state of E pro- 
vided that some (rank-1) projective measurements {11^} have 
been performed on S (whence "S'-discord" Vs) and have 
given the result j with probability pj. 

QD has recently attracted much attention due to its im- 
portance in distinguishing classical and quantum correlations 
II22I425II . and in particular in showing how a one-qubit quan- 
tum computer may work ll26ll . One can also find operational 
interpretations for QD, e.g., by exploiting the concept of quan- 
tum state merging luTll . Vanishing and time-invariant dis- 
cord in quantum states have been investigated in, respectively. 
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FIG. 1. (Color online) The initial state of the system (S) and envi- 
ronment (E) is a generally mixed state with vanishing discord. The 
state of the ancilla (A) is chosen such that it purifies gsE{0); i.e, 
QSEAiO) = \4>{0))sea{'4'{0)\- The state of the combined system 
(SEA) hence will remain pure in time. After state tomography to 
determine gsA{t), we are able to conclude whether the map £ (acting 
only on the system) is Markovian or non-Markovian. 



Refs. m and Ref. iH. 

Of particular importance to our goal in this paper, we re- 
mind that in Ref. OOll . it has been shown that, except for a 
"zero-measure" set of quantum states, vanishing quantum S- 
discord of the initial state of the system and the environment 
is a necessary and sufficient condition for the dynamics of the 
system to be CR 

Given a system-environment interaction Hssit), whereby 
the evolution U{ti,t2) = Texp[— i/^^ HsE{t)dt\, one can 

define as below a CP map £{t + S, t) on Qs{t) if(f) PsEJt) has 
vanishing 5-discord (up to the limitations stated in Ref. OOll ). 

gs{t + S) = Ti- EpsEit + S,t)QsE{t)UlE{t + S,t)] 

= £{t + S,t)gs{t), (3) 

for a given t,S > 0. One can constitute a concatenation of 
these maps as ft = £{T, T-6) . . . £{t+6, t) . . . £{S, 0), which 
by construction is an "infinitesimal divisible" CP map. Such 
maps have been proven to be always represented by a Marko- 
vian dynamics iSfl. It is evident that 11 = £{T,0). Thus in 
summary, we have the following result: 

Proposition 1: The dynamics of a system S interacting 
with an environment E is Markovian if(f) the instantaneous 
combined state of the system and the environment has vanish- 
ing QD, i.e., Vs [QsE{t)] = Vt > 0. 

This result, as is, provides a method which is hardly ex- 
perimentally measurable because it requires the knowledge of 
gsE{t) at all times, which is too demanding. To overcome 
this problem, we relax its strict condition and just focus on 
the sufficiency part. Specifically, we use a lower bound for 
QD which is independent of the information of the environ- 
ment. To this end, we use a monogamy relation for QD in a 
pure tripartite state qsea as 



Ef[gEA\ ='Ds[gEs] +Se\s^ 



(4) 



where Ef \gEA] is the entanglement of formation of the state 
of EA IBlll . and Se\s = S[gEs] ^ '^[e's] is the conditional 



entropy ||32ll . Since the total state gsEA is assumed to be pure, 
from the Schmidt decomposition we have S [gEs] = S [qa] , 
whence 

Ef [gEA] = Vs [gEs] +S[gA]-S [gs] . (5) 

Now let us purify the system-environment initial state 
Qse{0) by attaching an ancilla A which is not interacting with 
S and E. The total state gsEA (t) will remain pure in time be- 
cause SE evolves unitarily in time, while A does not evolve. 
In this case, noting that entanglement of formation is always 
nonnegative, Eq. (|5| yields 

VsigEsit)] >S[gsit)] ~S[gAit)] =: A^aW- (6) 

This relation implies that to recognize non-Markovianity of a 
quantum dynamics, it suffices to measure AsA{t)', an strictly 
positive value for Asa at some time i* is a signature of non- 
Markovianity of the dynamics at any later time. 

Proposition 1 can also be used to define a measure for non- 
Markovianity power of a quantum dynamics in an interval 
(0, r) as follows: 



^nm(t) 



1 



Vs[gsE{t)]dt. 



(7) 



A nonzero value for Pnm is a necessary and sufficient condi- 
tion for non-Markovianity of the associated quantum dynam- 
ics. Since again computing Vs [gsE{t)] is difficult, we pro- 
pose 



-Pnm(t) 



1 

27 



{\AsA{t)\ + AsA{t))dt, 



(8) 



as a lower bound for non-Markovianity power. 

Remark 1. It should be noted that our QD measure is 
not a trivial extension of the entanglement measure for non- 
Markovianity power of a quantum dynamics. One of the 
characteristics of QD is that it can increase under local op- 
eration and classical communication (LOCC) while entan- 
glement cannot. Separable states gsE = TliiPiSs ® Be 
(O < Pi < 1, J2iPi ^ l) ^^^ ^^ generated by LOCC from 
some initial product state, say, gsE{0) = |0)(0| ® |1)(1|. 
Since there exist separable states with nonvanishing QD 112011 . 
LOCC may increase QD. As an explicit example, consider the 
preparation of the system and ancilla in the vanishing QD state 
Qsa{0) = |(|0)(0| (E)Wi + |1)(1| (E) W2), where Wi and W2 
are two orthogonal density matrices. A Hadamard channel 

£{T)g=[l-p{T)]g+p{T)HgH, (9) 

with H = ((T^^' + ct('^^)/V2, acting on the system yields 



esAiT) = {£{T)<i^i)gsA{0) 



^{[i-p(t 



-pl+)(+l}®f«l 



-{[1-p{t)]\1){1\+p{t)H{-\}®W2,(10) 

in which |±) is the eigenvector of a^^\ Since 

\psa(t), gs(T)(g) 11] 7^ 0, gsA{T) is a state with nonvanishing 
QD in, while VsigsAiO)] = 0. 
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FIG. 2. (Color online) Concurrence 1I3TI1 (blue line) of qsa is de- 
creasing in every arbitrary time interval [i is in ft = 1 units]. Hence, 
deciding about Markovianity of the dynamics is not possible through 
this entanglement measure. On the other hand, Asa [Eq. {6}] (the 
red curve), is positive in t £ [0, 6], which shows that Ds [QSE{t)] is 
nonzero in this interval; hence, the process is non-Markovian. 



Remark 2. The positivity of the associated Choi- 
Jamilkowski state of a quant um p rocess has been used as a 
measure for its Markovianity fw\. To use this criterion, one 
needs to find £{t + e, t). But the problem is that in general 
£{t + e, t) is not a well-defined quantum process/map. When 
the evolution is non-Markovian, £{t + e, t) is not necessarily 
a map or identifiable by a physically conceivable tomography 
scheme. In fact, in the middle of the process at some time i», 
it may happen that gssit*) has nonzero quantum correlation. 
Hence, the evolution of the system cannot be interpreted as 
the application of a quantum map on gs- Note also that the 
number of independent parameters to be evaluated to identify 
the process between a t and t + e thoroughly is 0{d'*^), where 
d is the dimension of the Hilbert space of the system; while 
a complete set of quantum states has only 0{d^) independent 
members. This implies that process tomography for identifi- 
cation of £{t + e, t) is in general infeasible. In our approach, 
however, all one needs is just state tomography of gsE{t) in 
every point of time, which is well-defined (although of course 
difficult). Besides, if one use the lower bound proposed, our 
measure seems experimentally feasible. 

Remark 3. It should be clear that, although we use the result 
of Ref. IBOll . here we have set different goals and strategy. 
We deal with the problem of non-Markovianity and provide a 
measure thereof by using the concept of QD, while Ref. llsoll 
deals with a condition for CPT-ness of quantum maps. 



III. EXAMPLES 

In the following, we illustrate the main idea of this paper 
through two examples. 

Example I. Consider the Jaynes-Cummings model, in 
which the system of interest is a two-level atom interacting 
with a cavity mode through the Hamiltonian 




FIG. 3. (Color online) S-discord (red) of example II vs. time [in 
h = 1 units] for gsE{0) = |0){01 (p\0){0\ + (1 - p)|l){l|), 
where \(f>) = (|0) + \1))/V2 and p = 0.5. 5-discord has been 
obtained from Eq. (|5} in which \'ip(t))sEA = {UsE{t) ® ll)[|<;6) ® 
(y^|00)EA + yi — p|11)£;a)]. The blue curve shows concurrence 
between the system and ancilla vs. time. Since concurrence increases 
only after t* = 3, the non-Markovianity property of the dynamics 
for f < t« is not captured by using the entanglement measure of 
non-Markovianity. 



where <t+ = |1)(0| and a^ are, respectively, the raising 
operators of the system and environment. We assume the 
initial state of the atom (5) and cavity {E) to be a sepa- 
rable state esB(O) = gs(0) iS5 |a)_E(a|, in which |a) is 
a coherent state, and A = 1. For the case of £>s(0) = 
e|0)(0| + (1 - e)|l>(l|, the purification with a qubit an- 
cilla yields \i}}{Q))sAE = (V^|00) + VT^|11)) ® \a). 
When e = 0.2 and a = 5, the concurrence of the sys- 
tem and ancilla decreases in time. Hence, the entanglement 
measure of non-Markovianity does not imply a reliable clue 
about Markovianity/non-Markovianity of the dynamics. But 
since /S.sA{t) > 0, we are sure that the evolution is non- 
Markovian (see Fig. |2]i. Calculation of the lower bound on 
non-Markovianity power of the dynamics gives Pnm(IO) = 
0.094. 

Example II. Consider two qubits (system-environment) in- 
teracting through the Hamiltonian 



.(3) 



H=^{l+ar>+a\'>-ar<^','), 



A3) 



A3) (3). 



leading to the evolution operator 



in which (cr(o),cr(i) 






(12) 



(13) 



a(2),a(3)) 



H, 



SE 



A(cr (gjflf 



® a), 



(11) 



{t,X,Y,Z) are the iden- 
tity and the Pauli matrices. 

We show that if the initial state of the environment is diag- 
onal in the computational basis, the state of the environment 
does not change at any time, gsit) ~ ^£(0), and the sys- 
tem and the environment remain unentangled in time — see, 
e.g., Ref. PSil for a note on lack of back-action in revival 
of quantum correlations. One may infer that the evolution of 
the system should then be described by a quantum dynami- 
cal semigroup. However, we show that it is a non-Markovian 



transformation because a nonzero QD is produced during the 
evolution due to the interaction of the system and the envi- 
ronment. Assuming psb(O) — £'s(0) ® Qe^\ the Kraus 
operators of the dynamics of the system are Ao = ^/po 1 and 

A\ ~ sfPi^'^^" ^^' where pj = (j\qe\3)- To find whether 
the obtained time-dependent map £" (t) • =^aAj(t')- JVAt) 
is Markovian or not, one needs to calculate the ^-discord of 
QSE{t)- It is straightforward to show that 



QsE{t) = ^ Pj e 
ie{o,3} 



-»t^«'/2„„rm p^t'^'^'/a ^ |j)(j|,(14) 



'gs{0) 



which is a separable state. Given the environment initially 
in the state gE{0) = Po|0)(0| + pi|l)(l|, it is evident that 
0E{t) = te(0). 

Now we use a result of Ref. [231 that if a bipartite quantum 
state has zero S'-discord, then [qse, £'s ® l] = 0. Since in 
this example [gsE{t), Qs{t) ^ ^] 7^ unless [£is(0), cr'^'] — 
0, thus except for these special states, 5-discord is nonzero, 
and hence the dynamics of the system is non-Markovian (see 
Fig.E). 



IV. SUMMARY AND OUTLOOK 

An interesting question in open quantum system theory is 
that given the interaction Hamiltonian of a system and its en- 
vironment whether the induced dynamics on the system is 
Markovian or not. Several methods have been suggested to 
assess non-Markovianity of a dynamics. Here, based on the 
idea that underlying correlations play key role in a reduced 
dynamics, we employ quantum discord (QD) to define a mea- 
sure for non-Markovianity. To alleviate the need for the state 
of the environment in calculating QD in our scheme, we have 
considered an experimentally measurable lower bound for QD 
by using a monogamy property for tripartite quantum systems . 
This bound being nonvanishing provides a sufficient condition 
for non-Markovianity of a quantum dynamics. 



We also have discussed two examples to illustrate the re- 
sult of this study. In the first example, it has been shown that 
QD can signal non-Markovianity in a process which eludes 
the entanglement measure of non-Markovianity. The second 
example has featured an explicit quantum dynamics whose 
non-Markovianity is strictly due to the production of quantum 
correlation not in the form of entanglement. Here, the initial 
state of the system and environment was a product state and 
remained separable in time, while the state of the environment 
did not vary in time either 

Having powerful non-Markovianity measures at hand, one 
can think of various applications in quantum information the- 
ory. For example, an intriguing question is that how one 
should engineer an environment so that the generated sub- 
dynamics of the attached system becomes Markovian. This 
is important because, for example, it has been argued H that 
one can perform adiabatic quantum computation in an indirect 
fashion in an open system, provided that the effective system 
dynamics is a Markovian evolution with fixed points encod- 
ing ground states of some given Hamiltonian. Such scenar- 
ios require ability to decide whether an environment induces 
Markovian sub-dynamics. Our method provides a way to go 
to answer this question. Given an interaction Hamiltonian 
with unknown coupling constants. A, and a preparation of the 
system and environment (assuming being feasible), £>5_e(0), 
we can find the state of the system and the environment as a 
function of A, £>s_b(A, t), and check for which values of A, S'- 
discord remains zero in time. As a result, we are in principle 
able to manipulate coupling constants such that the desired 
map becomes Markovian or non-Markovian on demand. 
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